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Abstract 

We investigate a simple 5D extension of the Minimal Supersymmetric (SUSY) Standard 
Model (SM) that is combined with the bulk matter Randall-Sundrum (RS) model, which gives 
a natural explanation the Yukawa coupling hierarchy. In this model, matter and gauge super- 
fields reside in the 5D bulk while a SUSY breaking sector and the Higgs doublet superfields 
are localized on the infrared brane. The Yukawa coupling hierarchy in SM can be naturally 
explained through the wavefunction localization of the matter superfields. While sparticles 
obtain their flavor-blind soft SUSY breaking masses dominantly from the gaugino-mediated 
SUSY breaking, flavor-violating soft terms arise through the gravity-mediated SUSY breaking 
which are controlled by the wavefunction localization of the matter superfields. This structure 
of the model allows us to predict the sparticle mass spectrum including flavor-violating terms. 
We first explicitly determine the 5D disposition of matter superfields from the low energy ex- 
perimental data on SM fermion masses, CKM matrix and the neutrino oscillation parameters. 
Then, we calculate particle mass spectra and estimate the effects of the flavor-violating soft 
terms, which should be compared with the current experimental constraints. With gravitino 
being the lightest sparticle (LSP), the next-to-LSP, which is long-lived, is predicted most likely 
to be either singlet smuon-like or selectron-like. The model can be tested at collider experiments 
through flavor-violating processes involving sparticles. The flavor structure among sparticle, 
once observed, gives us a clue to deep understanding of the origin of Yukawa coupling hierarchy. 



1 Introduction 



The gauge hierarchy problem has been the main motivation for physics beyond the Standard 
Model. One notable solution to this problem is offered by Randall- Sundrum (RS) model PQ, 
which connects 4D Planck scale and the electroweak scale by means of 5D warped geometry. 
The RS setup also accommodates a natural explanation to the Yukawa coupling hierarchy in 
its extension with matters in the 5D bulk [2]. This bulk matter RS model solves the Yukawa 
hierarchy problem with the following common structure. We put Higgs on the infrared (IR) 
brane and fermions in the bulk. With non-hierarchical 5D Dirac masses, we localize the heavy 
fermions towards the IR brane and the light ones towards the ultraviolet (UV) brane. The 
wavefunction overlaps between the Higgs field and the fermion fields give rise to the Yukawa 
coupling hierarchy. In spite of their elegant solutions to the gauge hierarchy problem and the 
Yukawa hierarchy mystery, the RS models are under severe experimental constraints concerning 
the Kaluza-Klein (KK) scale due to the flavor- changing neutral currents caused by the KK 
states. Apart from the proton decay problem, the most stringent constraint comes from the 
data on K°—K° mixing and the 1st KK gluon mass should be > 21 TeV [3]. This constraint thus 
spoils their solution to the gauge hierarchy problem as well as their experimental accessibility 
at the Large Hadron Collider (LHC). (However imposing some flavor symmetry softens the KK 
scale bound, see [1].) 

In this paper, we study the supersymmterization of bulk matter RS model [5] . This model 
resorts to supersymmetry (SUSY) to solve the gauge hierarchy problem while maintaining the 
natural explanation of bulk matter RS model to the Yukawa hierarchy; the KK scale can be 
much higher than the electroweak scale and SUSY fills the gap between them. We especially 
investigate 5D Minimal SUSY Standard Model (MSSM) in RS spacetime, where we allow the 
gauge superfields to propagate in the bulk but localize Higgs superfields and the SUSY breaking 
sector on the IR brane, and the matter superfields are laid in the bulk with various 5D profiles. 

All phenomenological SUSY models require a realistic SUSY breaking mechanism in har- 
mony with experimental bounds. Our model naturally incorporates gaugino-mediated SUSY 
breaking mechanism as the matter superfields and SUSY breaking sector are (partly) sepa- 
rated by the 5th dimension while the gauge superfields couple to the SUSY breaking sector 
without suppression and mediate its effects to matter sector. At the same time, it explains 
the Yukawa coupling hierarchy through 5D localization of matter fields. Due to this structure, 
the sequestering is always incomplete; the superfields of 3rd generation particles lean towards 
the IR brane so that they couple to the SUSY breaking sector with less suppression. We thus 
have a unique pattern of flavor-violating soft mass terms that are related to the origin of the 
Yukawa hierarchy. Current experimental bounds on flavor violations give strong constraints on 
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the model and allow us to make predictions on the sparticle mass spectrum as well as on its 
flavor- violating effects. A similar setup was proposed in [6] as a 5D realization of "flavorful 
supersymmetry" [7]. 

Our model is a simple 5D extension of the MSSM in the RS spacetime, but has the ability 
to simultaneously provide a viable SUSY breaking mechanism and explain the hierarchical 
structure of Yukawa couplings. Since the sparticle mass spectrum and the Yukawa coupling 
hierarchy are rooted on the same 5D setup, this model has a strong predictive power on both 
of flavor- conserving and flavor-violating soft SUSY breaking terms. 

In the next section, we write down a general form of MSSM in the bulk of RS spacetime 
equipped with gaugino- mediated SUSY breaking. In section 3, we assume that all couplings in 
5D theory are of 0(1), and any hierarchical structure in 4D theory originates from 5D geometry. 
Based on this assumption, we determine 5D disposition of matter superfields from the observed 
fermion masses, Cabibbo-Kobayashi-Maskawa (CKM) matrix and neutrino oscillation data. 
In section 4, we make general remarks on the SUSY breaking mass spectrum. In section 5, 
we discuss the difference between our model and minimal flavor violation. In section 6, we 
calculate a sample of mass spectra and study experimental bounds on them. In section 7, 
we discuss a signature of the model, that is, unusual next-to-lightest sparticle (NLSP) and its 
flavor-violating decay. The last section is devoted for conclusion. 



2 Setup 

We consider 5D warped spacetime with the metric p]: 

ds 2 = e-^ri^dx^dx" - dy 2 , (1) 

where y is the 5th dimension compactified on the orbifold S 1 /Z 2 : —irR < y < ttR , and k is 
the AdS curvature that is of the same order as the 5D Planck scale M 5 . Assuming that the 
warp factor, e~ kRw , is much smaller than 1, we have the following relation for k and M5 : 

Ml = %l- e -^) ~ (2) 
k k 

where M* is the 4D reduced Planck mass. This relation implies k ~ M5 ~ M* . We put a UV 
brane at y = and an IR brane at y = irR. The fundamental scale on the UV brane is M5, 
while that on the IR brane is Mse~ kRn . Note that in our model, M§e~ kRT[ is not necessarily at 
TeV scale, but is at an intermediate scale between M* and TeV. 

All MSSM superfields reside in the bulk. However, for simplicity, the Higgs superfields are 
assumed to be localized on the IR brane. We adopt Polonyi model [S] for the SUSY breaking 



2 



sector and introduce a gauge-singlet chiral superfield X on the IR brane, whose F-component 
develops VEV to break supersymmetry. We take advantage of Giudice-Masiero mechanism |10j . 
namely, we impose an appropriate R-symmetry to forbid the SUSY-conserving /i-term and force 
the /i-term to arise from SUSY breaking effects. In this paper, we assign the following R-charges 
to X and the MSSM superfields: 

X : 0, H u/d : 0, QJU,/ D.J L,/ E t : +1, 

where H u , H d , Q iy Ui, D i: Li, E i: respectively, denote the chiral superfields of up- type Higgs 
doublet, down-type Higgs doublet, SU(2) doublet quark, singlet up-type quark, singlet down- 
type quark, doublet lepton, singlet charged lepton. Note that the above assignment permits 
higher dimensional superpotential for light neutrino Majorana masses. 

The 5D bulk action is described with 5D Af = 1 gauge multiplets and matter hyper- 
multiplets. We use 4D superfield formalism extended with the 5th dimension y, following [5]. 

An off-shell 5D Af = 1 gauge multiplet consists of a 5D gauge field A M (M = 0, 1, 2, 3, 5), 
two 4D Weyl spinors Ai,A2, a real scalar E, a real auxiliary field D and a complex auxiliary 
field F, all of which transform as the adjoint representation of some gauge group. They are 
composed into one 4D Af = 1 gauge superfield V and one 4D Af — 1 chiral superfield x t na t 
are 

v = -Qo^qa^ - ie 2 ex 1 + %d 2 d\ x + ^e 2 e 2 D , 

X = -^(Z + tA 5 ) + V29\ 2 + 9 2 F . 
V 2 

Under Z 2 parity: y — > —y, they transform as 

V -> V , x -> -x ■ 
The action for 5D Af = 1 gauge multiplets is given by 



S 5D „ = dy d 4 x e-^ 



l — J d 2 6e klyl trj (e^W aa )(e^W*) + h.c. } 



(9i) 



l — J d 4 0e 2fe Wtr{ (V2d y + x a ^ V (-V2d y + x a )e v - (d y e~ v )(d y e v )} 



(3) 



where a labels gauge groups and W aa denotes the field strength of V a in 4D flat spacetime. 
When the unitary gauge, A^ = 0, is chosen, only V a has a massless mode in 4D picture. This 
mode has no dependence on y and will be written as V (x, 9, 6). 
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A 5D M = 1 hypermultiplet is expressed in terms of two 4D M = 1 chiral superfields 
$, $ c that are in conjugate representations of some gauge group. We assume that the former is 
Z 2 -even and the latter Z 2 -odd. Taking the basis of diagonal bulk mass, we have the following 
action for 5D M = 1 hyper-multiplets: 



S 



5D chiral 



j dy j d 4 xe- 4%l J d 4 6e 2k ^ ($t e -^ + ^e v ^ h 



+ J d 2 6e kM Qftdy - X /V2 - (3/2 - a)k}$i + h.c. 



(4) 



where % is a flavor index and c$ denotes the 5D bulk mass in unit of AdS curvature k. Only $j 
has a massless mode in 4D picture, which will be written as ^)e^ 2 ~ c ^ fc ' s/ '. 

Let us write down the low-energy 4D effective action of the fields in the bulk, which is 
described with the massless modes of 5D M = 1 gauge multiplets and 5D M = 1 matter 
hyper-multiplets. After integrating over y, we obtain the following 4D effective action: 



S 



ADeff. 



d 4 x 



2ttR 



a , 2 



Jl~2ci)kR-K _ 

d 2 e w aa w a a + h.c. + / d 4 # 2— — 

(1 2c,- ) k 



>J e~ v c 



(5) 



where the dimensionful 5D gauge coupling, g§, is connected to 4D gauge coupling g\ by the 
relation: g§ = \/lTxRg%. 

Next we consider the theory on the IR brane. The IR scale, M^e~ hR7T , is a free parameter 
of the model and is only assumed at an intermediate scale between the 5D Planck and the 
electroweak scales. 

On the IR brane, we introduce Polonyi model for SUSY breaking: 



Sir D 



d A x 



d A 9 e 



-2kR.n 



d 2 e n\x + h.c. 



(6) 



where the term is for stabilizing the scalar potential of X at the origin, fix satisfies 

M 5 e~ kRn x TeV , (7) 



~2kRn,2 
e fl x ~ 



which is equivalent to 



M 5 



TeV 



M.e~ kR * ' 



so that it gives rise to gaugino masses at TeV scale through the VEV of Fx- Note that the scale 
of Hx is between the 5D Planck and the IR scales. This scale is put in by hand, as in tree-level 
SUSY breaking models, or is generated through a dynamical SUSY breaking mechanism [9] , of 
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which the Polonyi term (6) is the effective theory. We additionally assume that only the SUSY 
breaking term explicitly breaks the R-symmetry. 

Other terms on the IR brane are listed below: 

MSSM term: 



Sir D 



J d 4 x J d*6 e~ 2kR * { Hie~ v H u + E\e~ v E d } 
+ I d 2 9 e- 3kR * { e (3-*-«y)w&r ^k HUQ + e 

y i m 5 j 

J M 5 

Gaugino mass term: 



(3 — Cfc— C; )kRn 



(Vd) 



kl 



Mb 



HdDkQi 



Sir D 
Higgs SUSY breaking term: 
d 4 x 



d 2 fl + h.c. 



h.c. 

(9) 

(10) 



!fefl7r <! d mi ^r^H u H d + d bmu ——^-H u H d + h.c. 

M? 



+ / d 4 £ e~ 2fc ^ | d uA ^±^HlH u + du0 ^HtH u 



(11) 



Matter soft mass term 

ftfl D y d 4 a: 

+ (Q —y U, D, L, E ) . 
A-term-generating term: 



M| 



(12) 



/ d 4 x [ / d 2 8 e- 3kR * { e^-^ km ^IXHuUiQj + e <3-«*-«*)*«* i^L X H d D k Q 

y [y i M i M i 



We omitted brane kinetic terms because they only affect the overall normalization of the fields 
and are irrelevant to the point of our model. 

We normalize X, H u , H d , Qi, Ui, Di, Li, E\ to make their kinetic terms of the 4D effective 
theory canonical. This is done by the following rescaling: 



(13) 



X -»■ X = e~ kR7T X, H u -»■ H u 



-kRn 



H v , H d — y H d = e kR7T H, 



->■ 



e (l-2ci)kR-K _ l 

' (l-2a)k 



(14) 



where 0j denotes Qi, Ui, Di, Li or E^. Then the MSSM term becomes 



Sir D d 4 x 



y d 4 x y d*e { ift e -^^ M + #t e -^ | 



d 2 e 



1 - 2c- 



1 - 2c, 



fc 



+ 



1 - 2c k 



2{i - e -(i-2ci)fcj&r} y 2{1 - e -(i-acj)*ifer} m 5 

' k 



(jju)ij H u UiQj 



l-2c, 



2{1 - e -(l-2c fe )fcfe} V 2 {1 - e -(l-2 Ci )fcfe} jVf; 



(y d )fci H d D k Qi 



1 - 2c. 



l-2c n 



k 



2{i - e-a-^jfcflir} y 2 {1 - e -( 1 - 2c ") feR7r } Af 5 
The gaugino mass term will be 



(y e )mn H d E m L n > + h.c. 



S m D d 4 x 



d 2 # d, 



.(15) 



(16) 



The Higgs SUSY breaking term will be 



->IR 



j d 4 x 









M 5 e 



-kRw 



~ ~ x ] x ~ ~ 

H u H d + d bmu —K — -r^-H u H d + h.c. 



M 2 e -2kRn 



M 5 e- kR * u M|e- 2feR7r u 

+ ddA M^e~ kR ^ d ^ d + dm M 2 e~ 2kR ^ d ^ d 



The matter soft mass term will be 
SVr d / d 4 x 



(17) 



d 4 # 



1 - 2c, 



1 - 2c. 



M 5 e 



k 

2{i - e -(i-2 Cl )fci?7T} y 2 {i _ e -(i-2c,-)fcJ&F} W 5 



x 



+ (Q -4 C/, D, L, E 
The A-term-generating term will be 



Sir D / d 4 x 



d 2 9 



1 - 2c 



1 - 2c,- 



k (a u ) 



2{1 - e -( 1 - 2c i) fefe } y 2{1 - e "(i-2c,)fcfe} M5 M 5 e~ kRn 

k {ad)kl XH d D k Q t 

} V 2{1 - e-d" 2 ^)^} M 5 M 5 e- fc ^" XHdE ™ L * 



XHJJiQj 




h.c. 



(19) 
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We introduce light neutrino masses by simply writing down higher dimensional operators 
on the IR brane, namely, 



S IR D / d 4 x \ d 2 9 e~ 3kR * e ^- c ^ kR * (Y u ) pq LpHu ^ qHu + h.c. 
J J M 5 



1 2c p / 1 2c g / y \ LpH u L q H u 

a X I a V \l 2 {1 _ e -(l-2 Cp )fci?7r|y 2 {1 - e -(l-2c,)fcifcr} W« M^'^ 

(20) 

Note that the Kaluza-Klein (KK) scale, M^,e~ kRn , is related to the scale of light neutrino masses. 
Another possibility is to introduce singlet neutrino superfields and adopt the seesaw mechanism 
pTj . In this case, the KK scale can be a free parameter of the model. 
Now the MSSM Yukawa couplings are expressed as 



2{1 - e -(l-2c,)kRnj V 2 {1 _ e -(l-2 Cj )fci?7r j M , 
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l-2c fc 1 ~ 2c, fc_ 

A ' 2{1 - e^ 1-26 *)***} V 2{1 - e -(i-2ci)feR7rj.^ 5 



(^e)mn - \\ 2 {\ - e ~(l-2c m )kRnj \J 2 {\ - e ~(l~2c n )kRn} ^ (^)mn , (21) 

and the neutrino mass matrix m v is given by 



{m v ) pq - W 2 {i _ e -{iJc v )kR^ J 2 {i _ e -(i-i)w&r} ( Y »)p* M 5 e- kR ^ ' 



The geometrical factor a/(1 — 2c) / (2{1 — e"( 1 ~ 2c ) fci?7r }) has a unique property. For c < 1/2, it 
is approximated by — c and is 0(1). For c > 1/2, it is approximated by \/ c — 1/2 e -( c - 1 / 2 ) fcR7r 

and is exponentially suppressed. Therefore this factor can generate the large hierarchy of the 
Yukawa couplings without hierarchy. In the following, we assume that the components of 5D 
coupling matrices, y u , yd, y e , Y», are all 0(1) and that the hierarchical structure of MSSM 
Yukawa couplings and the neutrino mass matrix arises from the following terms: 



2c, ; / 1 - 2 Cj 



2{1 - e -{l-2a)kRK} V 2{1 _ e -(l-2cj)fcii7r\ 

We define geometrical factors ai as 



a; 



2 CT7^T with . = 1,2,3 (23) 



2{1 - e -(l-2c 9i )feB7r\ 
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for the i-th generation of SU(2) doublet quark superfields. Similarly, we define 7,, 5 iy 
for SU(2) singlet up-type quark, singlet down-type quark, doublet lepton, singlet neutrino and 
singlet charged lepton, respectively. Thus, the up- type quark Yukawa matrix Y u , the down- type 
Yukawa matrix Y d and the charged lepton Yukawa matrix Y e (in the basis of diagonal 5D bulk 
mass) are given by 

(Y u )ij ~ Pi<Xj , (Ydjij ~ liUj , (Y e )n ~ U&j , (24) 
and the neutrino mass matrix m v is given by 

v 2 

W« ~ ^ M 5 e- fc ^ ' (25) 
with VEV of the up-type Higgs doublet u u . 



3 Yukawa coupling hierarchy from geometry 

In this section, we determine the order of the geometrical factors, a», Pi, 7$, 5$, 6j, from the 
experimental data on SM fermion masses, CKM matrix and the neutrino oscillation parameters. 
Note that the geometrical factors must be evaluated at the KK scale, ke~ kRn , where the 5D 
theory is connected to the 4D effective theory. However, as is seen from [12], the renormalization 
group (RG) running changes the Yukawa couplings by at most a factor 2 and CKM matrix 
components by at most 1.2 through the RG running from ~ 10 15 GeV to electroweak scale. 
Also the neutrino mass matrix is affected only by 0(1) through the RG running [13]. Therefore 
we can estimate the order of a*, Pi, 7$, 5i, e, directly from the experimental data at low energies. 

We first show the model's predictions on Yukawa eigenvalues and CKM matrix. Let us 
diagonalize the Yukawa matrices: 

V u Y u Ul = diag , 
V d Y d U\ = diag, 
V e Y e U\ = diag. 

For successful diagonalization of the hierarchical Yukawa matrices, the unitary matrices, 
U u , Ud, V u , Vd, U e , V e , need to have the following structure: 

/I 0\ 

U u ~ U d ~ a 1 /a 2 1 , V u ~ (a-tP), V d ~ (a -> 7) , 
\ ai/a 3 a 2 /a 3 1 / 

U e ~ (a -»• 5) , K ~ (a ->• e) (26) 
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which leads to 



V u Y u lll ~ diag ( ftai, (3 2 a 2 , /3 3 a 3 ) , 
V d Y d U\ ~ diag ( 71 «i, 7 2 a 2 , 73«3 ) , 

KW] ~ diag ( eA, e 2 5 2 , e 3 5 3 ) • (27) 

The hierarchical structure of CKM matrix Uckm is given by 

(1 cti/a 2 a>i/o>z \ 
«i/a 2 1 ct 2 /a 3 . (28) 
«i/a 3 a 2 /a 3 1 / 

The absolute values of the CKM matrix components, \Uckm\, at electroweak scale has been 
measured as PET 



\U C km[M w ]\ 

0.97419 ±0.00022 0.2257 ± 0.0010 0.00359 ± 0.00016 

0.2256 ±0.0010 0.97334 ± 0.00023 0.0415 ± 0.0010 - 0.0011 

0.00874 ± 0.00026 - 0.00037 0.0407 ± 0.0010 0.999133 ± 0.000044 - 0.000043 

We approximate this matrix by the following formula: 

/I \ A 3 \ 

\U C km\ ^ A 1 A 2 with A = 0.22 . (29) 
V A 3 A 2 1 J 

To discuss the neutrino mass matrix, we adopt the tri-bi-maximal mixing matrix [15] (which 
gives almost the best fit in the neutrino oscillation data): 

\ 

3 Y 3 u 

1 fi PL 

6 V 3 V 2 

/i R fi 

\ ~Ve V 3 ~V 2 / 
and the following data on neutrino mass squared differences [H] : 

Am 2 21 = 7.59 ± 0.20 x 10" 5 eV 2 , |Am 2 2 | = 2.43 ± 0.13 x 10~ 3 eV 2 . 

Also we assume that the mass of the lightest neutrino is negligible, for simplicity. Then the 
neutrino mass matrix, Umns diag(m„i, m v2 , m vz)UM N s ) is given by 

0.29 0.29 0.29 

U MNS dmg(m ul ,m u2l m u3 )Ul INS = ( 0.29 2.8 -2.2 ] x KT 11 GeV (30) 

0.29 -2.2 2.8 



U 



MNS 



for the normal hierarchy case, while 

/ 4.9 0.026 0.026 \ 

Umns dmg{m vll m v2 ,m^)Ul INS = 0.026 2.5 2.5 x 10^ n GeV (31) 

\ 0.026 2.5 2.5 / 

for the inverted hierarchy case. 

Now we are ready to compare the model parameters with the experimental data and estimate 
the order of oti, fy, ji, Si, 6j. For Yukawa eigenvalues, we simply have 

P\oti ~ m u /vsinf3 , P2&2 ~ m c /vsmf3, faces ~ m t /vsmf3, (32) 
7i«i ~ md/vcos(3, 7 2 a 2 ~ m s /vcos(3, 7303 ~ rrib/vcos(3, (33) 
€i5\ ~ m e /vcos(3, €262 ~ m^/vcos(3, e 3 5 3 ~ m T / v cos f3. (34) 

Since the top Yukawa coupling is of ~ 1, we have a 3 /3 3 ~ 1, which leads to 

«3 ~ 03 ~ 1 • (35) 

Comparing (28) with (29), we find 

«i ~ A 3 , «2 ~ A 2 . (36) 

We then have 

/?i ~ A~ 3 m u /vsm(3 , /3 2 ~ A~ 2 m c /t;sin/3 , (37) 
71 ~ A~ 3 m d /v cos /3 , 72 ~ A -2 m s /vcos[3 , 73 ~ m b /vcos(3. (38) 

Next compare the matrix (25) with the observed neutrino mass matrix. For the normal 
hierarchy case, it is possible to reproduce the hierarchical structure of the neutrino mass matrix 
by adjusting 

35± ~ $2 ~ $3 (39) 

with the factor 3 coupling of the 5D theory. On the other hand, for the inverted hierarchy case, 
we cannot reproduce the neutrino mass matrix with 0(1) couplings. The situation gets worse 
if we consider non-negligible mass of the lightest neutrino. Therefore, the model favors the 
normal hierarchy of neutrino masses with the relation (39). We estimate q from the relation 
(39) as 

61 ~ 3 5$ m e /v cos/3 , e 2 ~ S^ 1 m fl /vcosf3 , e 3 ~ S^ 1 m T /v cos/3 . (40) 
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Finally, we refer to the connection between the light neutrino mass scale and the KK scale. 
If the neutrino mass arises from higher dimensional superpotential, as in (25), the two scales 
are related through the following formula: 

3 x KT 11 GeV . (41) 



3 M 5 e- fci ^ 

Based on the relation above, we can estimate the KK scale from the value of 5%. 



4 Two origins of soft SUSY breaking terms 

In this model, SUSY breaking terms have two origins. One is contact terms between the SUSY 
breaking sector and the MSSM sector on the IR brane (gravity mediation contributions) |16j . 
The other is radiative corrections, in particular, the renormalization group effects from gaugino 
soft masses (gaugino mediation contributions) [T7]. For the superpartners of matter particles, 
the former induce flavor-violating soft terms while the latter mainly generate flavor-diagonal 
terms Due to the model's structure, the gravity mediation contributions are related to the 5D 
disposition of matter superfields that gives rise to the Yukawa coupling hierarchy. 

As is argued in [18], when the square root of space-like momentum, p = a/— p 2 , is larger 
than the KK scale, ke~ kRlT , 5D gaugino propagator connecting the UV and the IR branes is 
suppressed by the factor 

exp[ -p/{ke~ kR *) ] . 

Since the integral of loop momentum is done with Euclidized momentum, l 2 E = —I 2 , loop dia- 
grams containing gaugino propagators between the two branes are also exponentially suppressed 
when the range of integral is limited to [0(ke~ kR7T ), 00), that is, when the renormalization scale 
is around the KK scale. Matter superfields confined on the UV brane receive SUSY breaking 
effects through loop diagrams involving gaugino propagators in the bulk and gaugino mass on 
the IR. Hence we argue that, at the scale of ke~ hR7T , matter SUSY particles in the bulk gain 
SUSY breaking mass only through the contact terms on the IR brane, and radiative corrections 
through gauginos are negligible. When p < ke~ kRn , 5D gaugino propagator approaches to the 
4D one divided by irR, and matter SUSY particles receive SUSY breaking effects through gaug- 
ino radiative corrections just as in 4D MSSM. Based on the discussions above, we calculate the 
SUSY breaking mass spectrum in the following way: At the renormalization scale fi r = ke~ hRn , 
SUSY breaking terms arise only from contact terms on the IR brane (gravity mediation). In 
particular, 1st generation matter sparticles that are localized towards the UV brane have al- 
most zero soft mass. Below the scale of ke~ hRn , the RG equations of 4D MSSM controls the 
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mass spectrum (gaugino mediation). Therefore we can calculate the sparticle mass spectrum 
at the electroweak scale by solving the MSSM RG equations with the initial condition that, at 
\i r = ke~ kRn , SUSY breaking terms be given by the IR brane contact terms. In the rest of the 
paper, we denote the scale ke~ kR7T as M cut . 

At \i r = M cut , the SUSY breaking terms are given as follows: 



gaugino masses My 2 
Higgs B/i term B^i 
Higgs soft masses 



- d a A{gtf <Fx> 



dbmu 

{—duo 



M 5 e- 

< Ex > I 2 



kRir 



d 2 uA ) 



-d d0 + d 2 dA ) 



\<F X >? 
M 2 e~ 2kR * 

< F x>\ 2 



Mle- 2kR « 



matter soft masses 



m 



QJij 



-dQ2ij + d QUj ) onay 



<Fx>\' 
M 2 e~ 2kR * 



(Q, a) 
A — terms (AJy 



(U,/3), (D, 7 ), (L,5), (E,e) 

( \ _ <F x > 

-d uA (Vuhj PiVj — M5e _ kRn 



[a 



(42) 
(43) 
(44) 
(45) 

(46) 

_k_ < Fx > 

3 M 5 M 5 e~ kR « 



<F X > 
HJ M 5 e- kR * 

-ddA (Ydjij 



+ {a u )ij fra 



k 



<F X > 



<F X > 



-d d A (Y e ) 



M 5 e~ kRn 
<F X > 
M. e - kRn 



+ (a>d)ij JiOij 



1 M 5 M 5 e~ kR * 
k < F x > 



M 5 M 5 e~ kRn 
k < F x > 



(47) 
(48) 
(49) 



where ai, Pi,ji,5i,ei are defined as in (23). In addition, the //-term arises from the SUSY 
breaking effects (Giudice-Masiero mechanism): 



d r , 



<F X > 
M,e- kRn 



(50) 



Note that the flavor structure of matter soft masses and A-terms corresponds to the Yukawa 
coupling and neutrino mass matrix hierarchy in a unique way, governed by Pi, ji, Si, e^. 

We solve the MSSM RG equations from M cut toward low energies with the initial conditions 
(42-50), and evaluate the sparticle mass spectrum at the electroweak scale. 

Finally we remark on the nature of the lightest SUSY particle (LSP) and the next-to-lightest 
SUSY particle (NLSP) in this model. The gravitino mass is given by 



m 



3/2 



<F X > 



<F X >\ M 5 e 



-kR.Tr 



~kR.iT 



TeV x e 



-kRm 



(51) 
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and thus gravitino is always LSP, as in [19]. NLSP mainly consists of singlet sleptons whose 
flavor composition depends on the amount of gravity mediation contributions. Normally, the 
singlet stau is lighter than smuon and selectron due to its large Yukawa coupling, but in this 
model, it gains large soft mass through gravity mediation and may not be the lightest. As 
with other gravitino LSP scenarios, NLSP is long-lived because its coupling to gravitino is 
suppressed by 1/| < F% > \. 

5 Comparison with Minimal Flavor Violation 

The minimal flavor violation (MFV) is the setup that only SM Yukawa couplings violate flavor 
symmetry. In MFV, flavor-violating soft terms are generated via the MSSM RG equations 
involving Yukawa couplings. We here estimate the orders of the flavor-violating soft terms 
generated through RG running in MFV, and compare them with those via the gravity mediation 
in our model. We will see that the latter show different patterns from the former. 

We first introduce a flavor basis where Y u or Yd and Y e are diagonalized by the following 
unitary matrices £/*: 

UijY u Uq u = (diag.) , 
U D YdU Qd = (diag.) , 
U E Y e U L = (diag.) . 

Note that U*s depend on the renormalization scale as Yukawa matrices receive RG corrections. 
We will estimate the orders of the changes of f/*'s through RG running. Below is the RG 
equations for Y u : 

^(UuY u U Q u) = (^Uu^UuYMqu) + Uu(fi-^Y u )U Qu + (UuY u U Qu )U j Qu (fi-^U Qu ) 

= (^Uu^UuYJJqu) 
d/x 

+ 16^2 Uu { YuY d Y d + WuY*Y u + 3tr[YjTJ Y u + ti[Y ] D Y D ]Y u 

-(^9l + 3g! + fgl)Y u }U Qu 
+ (UuY u U Qu )U^ Qu (^U Qu ) , (52) 

where Yjj is neutrino Dirac Yukawa coupling which appears when we introduce the singlet 
neutrinos lighter than M cut . We hereafter adopt the GUT normalization for g±. From (52), we 
see that U\jY u Uq u remains diagonal during RG running when the unitary matrices satisfy the 



13 



following conditions, 

H^-Uu = 0, (53) 
an 

V-^Uq u = -^^2 ( off ~~ diagonal components of Y^Y d ) U Qu . (54) 

In the same manner, we obtain the following conditions for keeping UoYdUQd and UeYJJl 
diagonal: 

p—U D = 0, (55) 
an 

fJ,—U Qd = T"7^~ 9 ( off ~ diagonal components of Y^Y U ) U Qd , (56) 

d/i l07T' s 

li— U E = 0, (57) 
an 

u-^-Ul = t; (off — diagonal components of YXY D ) XJ l . (58) 

d/i 16n 2 

Now that we know how ^-diagonal basis, Y^-diagonal basis and y e -diagonal basis change 
through RG running, we estimate the orders of MFV effects on A-terms in these bases. The 
MSSM RG equations for A-terms are given by: 

lQn 2 n^-A u = 3A U Y^Y U + 3Y U Y^A U 



d/i 



+ A u YjY d + 2Y u YjA d 

+ 2( 3trptf4J - ^fiMff - 3g 2 2 M?J? - y^K/f )Y U 

+ ( 3tr[FjF u ] - ^g 2 _ 3 gl - fg 2 3 )A U 

+ tv[YiY D ]A u + tv[Y} ) A D }Y u , (59) 



lQ n ^± Ad = 3A d YjY d + 3Y d Y}A d 



+ A d Y^Y u + 2Y d Y^A u 

+ 2( 3tr[yjA d ] + tr[y e U e ] - ^M^ 1 - 3g 2 2 M?J? - y^K/f 

7 1 6 

+ ( 3tr[F d V d ] + tv[YX] ~ ~ ~ ysg )4, , (60) 



167r 2 /i-^-A e = 3A e y e f F e + 3Y e Y ] e A e 



d/i 



+ 2( 3ti[YjA d ] + tr[y e U e ] - ^<ftM$ - 3g 2 M?J 2 2 )Y e 
+ ( 3tr[yjy,] + tr[YjY e ] ~ \g\ ~ )A 

+ A e y^y D + 2Y e Y ] D A D , (61) 
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where neutrino Dirac Yukawa coupling Yd and the corresponding A-term Ad appear when we 
introduce singlet neutrinos lighter than M cut . From (53-61), we obtain the following equations 
for A u , Ad, A e respectively in Y u , Y d , Y e -bases: 

lQ7r 2 ^(UuA u U Qu ) = 3U u A u YJY u Uqu + 3U u Y u YJA u Uqu 



d/j, 



+ (UuA u U Qu ) (diagonal part of tf Qu Y}Y d U Qu ) + 2U u Y u Y}A d U Qu 
+ 2( 3tr [FjAJ - ^g£M$ - 3g\M{j 2 2 - ™g*M$ )(UuY u U Qu ) 

+ ( 3tr[Y^Y u ] - - 3g\ - ^gj ){UuA u U Qu ) 
15 3 

+ tT^YdiUuKUqJ + tAyhMiUuYuUQv) , (62) 



167T 2 fi^-(U D A d U Qd ) = 3U D A d YjY d U Qd + 3U D Y d YjA d U Qd 



dfi 



+ (U D A d U Qd ) (diagonal part of U ] Qd Y^Y u U Qd ) + 2U D Y d Y^A u U Qd 

7 16 
+ 2( 3tr[YjA d ] + tr[Y e U e ] - -gfM^ 1 - 3g 2 2 M^ 2 - -g 2 3 M?J? ){U D Y d U Qd ) 

7 16 

+ ( 3tr[Yjy d ] + tr[Y e tY e ] - -g\ - 3g\ - -g\ ){U D A d U Qd ) , (63) 



167T 2 ^(U E A e U L ) = 3U E A e Y}Y e U L + 3U E Y e Y}A e U L 
d/i 

+ 2( 3tx[YlA d ] + tr[Y e U e ] - j^M^ 1 - 3g 2 2 M^ 2 ){U E Y e U L ) 
+ ( 3tr [yJy,] + tr[YX] ~ \gl ~ Zgl )(U E A e U L ) 

+ (U E A e U L ) (diagonal part of U^Y^YdUl) + 2%F e F^ D f/ i , (64) 

To study the effects of MFV, we set the initial conditions for A u , A d , A e as 

{A u )ij lini. = M u (Y u )ij , [A d )ij | ini. — M d (Y d )ij , (A e )ij\i n i = M e (Y e )ij 

with mass parameters, M u , M d and M e . Then the terms 2UjjY u Yj A d XJq u in (62), 2U DY d Y^AJjQ d 
in (63) and 2U E Y e Y D A D U l in (64) respectively give rise to off-diagonal terms of (UuA u Uq u ), 
(UDA d UQ d ), {U E A e UL), which were initially diagonal. These off-diagonal terms in turn gen- 
erate off-diagonal terms through (62-64), but this does not change the orders of themselves. 
Noting that the orders of the Yukawa components in each basis are given as (Sy is the ordinary 
Kronecker's delta): 

(UuY u UQ U )ij ~ fiiOii Sij , (UuY d UQ U )ij ~ 7i«j , 
{U D Y u UQ d )ij ~ fiioij , (U D Y d UQ d )ij ~ -fiOti , 
(U E Y e U L )ij ~ , (U E Y D U L )ij ~ QSj , 
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where (j's are the geometrical factors for singlet neutrinos satisfying Q < 1, we estimate the 
orders of the off-diagonal terms of (UuA u Uq u ), (UdAJJqJ), (UeAJJ l ) that arise through RG 
running as (i ^ j): 

AiUuAuUQu)^ 



&{U D A d U Qd )ij 



A(U E A e U L )ij 



where M seesaw indicates the mass scale of the singlet neutrinos if they exist. We used the 
approximation that Efc(7fc) 2 - (73) 2 , Efc<A) 2 - (AO 2 and Efc(Cfe) 2 - (Cs) 2 - As for diagonal 
terms of (UjjA u Uq u ), (UdAJJqJ), (UeAJJl), the equations (65-67) do not change their orders. 
In conclusion, the orders of MFV effects on A-terms are given by the estimates (65-67). 

We next estimate the orders of MFV effects on matter soft mass terms, m 2 Q, mfj, m 2 D , m L , 
m 2 E , in the basis where Y u or Y d and Y e are diagonal. Below is the list of those terms in MSSM 



167T 2 

2 

16vr 2 
1 

16vr 2 
2 

16vr 2 
1 

16V 2 
2 

16tt 2 



In 
In 
In 
In 
In 
In 



M 



cut - 



2{UuY u YjA d U Qu ] 



ij 



Mi 

— ^) A«i a^faj M u , 

M w 

M cut , 



M 



w 



2{U D Y d Y^A u U Qd )i 



ij 



—^) jiai ai((3 3 ) 2 aj M d 
M cut 



M,, 



2(U E Y e Y^A D U L ) 



M. 



(65) 



(66) 



(67) 
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RG equations that give rise to flavor non-universal soft masses: 



\§ix 2 \i—m 2 Q 
d/i w 


D 


Y^Y u m 2 Q + m 2 Q YlY u + 2^m^y„ + 2(Y^Y u )m 2 Hu 




+ 
+ 


Y\Y d m 2 + mjyjy d + 2Y ] d m 2 D Y d + 2{Y ] d Y d )m 2 
24 A„ + 2A\A d , 


(68) 


2 d 2 
16n u—rriTj 
d/i 




2Y u Yjm 2 u + 2m 2 u Y u Y} l + AY u m 2 Q Y^ + 4(Y u yJ)r< 




+ 


4-A u 4 ' 


(69) 


2 d 2 
167r n——m D 
d/i 


D 


21^7^ + 2m^y d yj + AY d m 2 Q Yj + A(Y d Y})m 2 Hd 




+ 




(70) 


2 d 2 

1 (\lY II Tfl r 

d/i 


/ 


Y^Yml -\-m 2 -Y J( Y 4- 2Y^mlY + 2(Y i <Y)m 2 -, 

1 e 1 e< IL L * IIL L 1 e e ^ e " b E 1 e T" ^-f e 1 e)l IL H d 




+ 


2A\A e 




r d/i E 


+ 


Y D Y D m 2 L + mi^F D + 21^^ + 2(^F D )m^ + 2^+,^ , 


(71) 


D 


2Y e y e t m| + 2m|yy e t + 4^™^ + 4(y e y e t )m^ d 




+ 


4A e 4 , 


(72) 



where again, Yjj and Ad appear when singlet neutrinos lighter than M cut exist. From (53-58) 
and (68-72), we obtain the following equations for rnq in F u -diagonal basis, in Y^-diagonal 



basis, m 2 D in Y^-diagonal basis, m 2 L in Y" e -diagonal basis and m 2 E in Y" e -diagonal basis: 
! d ( 



16^ d ( , ^ uX u Y^Y u mQUQ U + UQ U mnY^Y u UQ U 



+ 2UlXm 2 v Y u U Qu + 2(U^ Q XYuU Qu )m 2 Hu 

+ ^(diagonal Parts of Y\Y d )UQ U (U ] Qu m 2 Q U Qu ) 

+ ( U Qu m Q U Qu)U f Qu (dmgon&\ parts of Y ] d Y d )U Qu 

+ 2U^mlY d U Qu + 2(U^}Y d U Qu )m 2 Hd 

+ 2U^ Qu AiA u U Qu + 2U ] Qu A\A d U Qui (73) 



l^H^-iUumlul) D 2U u Y u Y^m 2 u U{ J + 2U u m 2 u Y u Yjul 



dfx 



+ 4UuY u m 2 Q Yjul + A{UuY u Y^Ul)m 2 Hu 

+ AU v A u Aiul, (74) 



16n 2 ^(U D m 2 D Ul) D 2U D Y d Y ] d m 2 D U ] D + 2U D m 2 D Y d Y ] d Ul 



+ ttJ D Y d m%Ylu], + ±{U D Y d Yiu),)m 
+ 4U D A d A d Ul, (75) 



" l H d 



lQ7r 2 fi^-(Ulm 2 L U L ) D UlY}Y e m 2 L U L + Ulm 2 L Y}Y e U_ 



d/x 



L 

+ 2Ui^m%YJJ L + 2{UlY}Y e U L )m 2 Hd 
+ 2UlAlJU7 L 

+ C/J (diagonal parts of YlY D )U L (Ulm 2 L U L ) 
+ (t/lm|C/L)C/i(diagonal parts of Yj^)^ 
+ 2^yJm^C/ L + 2(C/i^y D C7 L )m^ + 2UlA* D A D U L , (76) 



167r 2 /^(£4;rr4f4) D 2U E Y e Y^m 2 E U E + 2U E m 2 E Y e Y^U E 



dfi 



2 



+ W E Y e m 2 L Y}ul + A(U E Y e Y^U E )m 
+ 4f/ E A e 44. (77) 

To extract the effects of MFV on matter soft mass terms, we consider the case with the following 
initial conditions: 



| ini. 


= ™*o% (* 


= Q,U,D,L,E,N) , 




{.Au)ij | ini. 


= M u (Y u ) l3 , 


{A d )ij ini. = M (J (F d )jj , 


(^4e)j;|ini. -^eO'eJij 


{A DJij | ini. 


= M D {Y D )ii . 







In (71), the terms involving or A d induce the off-diagonal terms of (U^itiqUqu) of the order 
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{i * h 



^ Qu m 2 Q U Qu ) 



,M 



ln(^) a^sfaj { 2m 2 + 2m D0 + 2m 2 Hd + 2M 2 } . (78) 



16tt 2 '" y M w 



l Q0 



It is then clear that, in Y^-diagonal basis, (^Qd m Q^<2 rf ) obtain the off-diagonal terms of the 
order (i ^ j): 



16tt 2 



ln( 



M 



cut ■ 



Mi 



w 



OiWaj { 2m 2 QQ + 2m 2 U0 + 2m 2 Hu + 2M 2 } . (79) 



If singlet neutrinos lighter than M cut exist, the terms in (74) involving Yd or Ad induce the 
off-diagonal terms of {Jj\m\]Ji) of the order {% ^ j): 



A(Ulm 2 L U L ] 



1 - , M cut 



167T 2 



ln( 



M,, 



■) S^)^ { 2m 2 L0 + 2m 2 m + 2m 2 Hu + 2M 2 D } . (80) 



On the other hand, for (UumfjU^), off-diagonal terms arise from the following two terms in 
(74): 

AUuY u m 2 Q YM = 4(UuY u UlJ(U Qu m 2 Q UlJ(UlXul) , 

WuAJivl = 4(Mu4)(t/L4^) ■ (81) 



Off-diagonal terms of (^q^q^qw) induced through (73) and those of (UuA u Uq u ) induced 
through (62) in turn give rise to off-diagonal terms of (UumfjUjj) through (81). Therefore, 
from (78) and (65), we estimate off-diagonal terms of (UumfjUjj) induced by RG running as 



2 rrt- 



UJ*3 



1 H M t . 



16tt 2 v M w 

2„, „, a /o™2 i o™2 i o™2 



x [ 4/3 l a l a J ( 73 )'a i a J /3 i (2m^ + 2m D0 + 2m 2 Hd + 2M 2 ) 
+ 4^ ( (3iai5 ik M u + ( 73 ) 2 a fc M u )( a k (3 k 5 kj M u + a fe (7 3 ) 2 a j a i /3 i ) ] 



1 . M cu t- 



ln(- 



167T 2 V Mi 



w 



x [ 4f3 l (a i ) 2 ( l3 ) 2 (a J ) 2 f3 j (2mJ + 2m^ + 2m 2 Hd + 2M d 2 ) 

+ 8A(a i ) 2 (7 3 ) 2 («,) 2 /3 J M 2 + 4A(«,) 2 (73) 2 («3) 2 (73) 2 («,) 2 ^M 2 ] . 



52) 



Likewise, we obtain the following estimates on off-diagonal terms of (UDfn 2 D U D ) and (L^m^L^ 



,2 rrt 
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{i * h 



A(U D m 2 D Ul 



ln( 



M, 



cut • 



A(U E m 2 E Ul 



X 

+ 

X 

+ 



_16tt 2 "~"M w 
4%(« J ) 2 (/33) 2 (a,)S (2mJ 

\1l a \2 1 „ \2_, 7i jr2 



2m 



uo 

2/ 



2m, 



\2f„ \2, 



2Mi 

\2i 



7i (« J ) 2 (/3 3 ) 2 («,) 2 7,Mi + ^(aOWM W(«i) V« ] , 
M cut 



53) 



ln(- 



167T 2 v M seesatu 

46 J (5 i ) 2 (C 3 ) 2 (^) 2 e, (2mi + 2m 2 V0 
Je*(*) 2 (Cs) 2 (*i)%iW3 + 46 4 (5 4 ) 2 (C3) 2 (53) 2 (C3) 2 (^) 2 ^M ( 



2m 



2M 



54) 



In summary, the orders of MFV effects on soft mass terms are given in (79, 80), (82-84). 

We have estimated the orders of MFV effects on A-terms and soft mass terms in the basis 
where Y u or and Y e are diagonal. In the rest of the section, we compare MFV effects with 
flavor- violating gravity mediation effects of our model and discuss their difference. 

Flavor- violating gravity mediation effects at the scale M cut can be read from (46-49). We 
assume that the couplings in 5D theory, d*, a*, are all 0(1). This is a natural assumption 
because we are trying to explain the hierarchy of 4D theory from 5D geometrical point of view. 
In an arbitrary basis, the flavor-violating parts of A-terms that arise from gravity mediation 
and are not proportional to the corresponding Yukawa couplings are given by 



{A d )ij D 
(A e )ij D 

at the scale M cut , where Mx is defined as 

M X EE 



eiSjMx 



I <F X > | 
M. e - kRjT 



(85) 
(86) 
(87) 



Matter soft mass terms that arise from gravity mediation are given by 



aiC(jM x , 



(Q, a) -+ (U,/3), (D, 7 ), (M), (E, e) (88) 

at the scale M cut . First, we argue that, at the electroweak scale, the flavor-violating parts of 
A-terms that are not proportional to the corresponding Yukawa couplings are still estimated 
as in (85-87) and flavor- violating parts of matter soft mass terms are estimated as in (88). 
This is understood from the form of RG equations; the right hand sides of the RG equations 
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(59-61) depend on the Yukawa couplings and A-terms themselves. For the component (A u )ij, 
the right hand side of (59) is at least proportional to fliOij. Similarly, the right hand sides 
of (68-72), which express flavor- violating contributions, depend on the Yukawa couplings and 
A-terms. The flavor- violating part of the RG equation for the component (mg)jj is at least 
proportional to ajjOj. The same discussion applies to other A-terms and matter soft masses, 
and we conclude that RG running keeps the orders of flavor-violating parts of A-terms and 
matter soft mass terms as in (85-88). Second, we argue that the estimates (85-88) are valid 
even in Y u or and Y e -diagonal basis. This is because the 5D couplings (y*)ij, d*uj and 

d*2ij in (46-49) are independent of each other. Therefore the matrices and d^ij, 

which give rise to A-terms and matter soft masses, are arbitrary even when (y*)ij is diagonal. 

Let us compare the orders of MFV effects (65-67, 79, 80, 82-84) and those of flavor-violating 
gravity mediation effects (85-88). We assume that M u , Md, M e in (65-67) and ^ 2 o> m jM 
in (79, 80, 82-84) are of the same order as Mx in (85-88). For the A-term components (A u )ij, 
(A u )2j in Y u -diagonal basis and (Ad)ij, {Ad)2j in ^-diagonal basis, MFV effects are always 
much smaller than flavor-violating gravity mediation effects because the former are suppressed 
by (ai) 2 or (c^) 2 compared to the latter. For the components (A u ) 3 j, (A^jsj, MFV effects can 
be of the same order as flavor- violating gravity mediation effects. For the A-term (A e )ij, MFV 
effects are much smaller than flavor-violating gravity mediation effects when the order of S3 is 
significantly smaller than 1 (we will see in the next section that this is the case for a realistic 
mass spectrum). If singlet neutrinos lighter than M cut do not exist, A e is diagonal. 

For the soft mass terms m 2 Q and m 2 L (if singlet neutrinos lighter than M cut exist), MFV 
effects can be of the same order as flavor- violating gravity mediation effects. For the terms mfj 
in Y u -diagonal basis and m 2 D in Y^-diagonal basis, MFV effects are always much smaller than 
flavor- violating gravity mediation effects except their (3,3)-components. This is because the 
components of these terms other than (3,3) are at least suppressed by («i) 2 or (o^) 2 - For the 
components (mfj)^ in Y u -diagonal basis and (7^0)33 in Y^-diagonal basis, the former can be as 
large as the latter. For the term m 2 E in Y e -diagonal basis, MFV effects are much smaller than 
flavor- violating gravity mediation effects when the order of £3 is significantly smaller than 1. If 
singlet neutrinos lighter than M cut do not exist, there is no MFV on m 2 , m 2 E . 

In this section, we have discussed the difference between the flavor-violating soft terms in 
MFV and those generated by the gravity mediation of our model. We have proved that, for some 
components of A-terms and soft mass terms, the gravity mediation contribution dominates. 
Therefore, it is in principle possible to distinguish our model from other SUSY models with 
MFV. 
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6 Particle mass spectra and experimental constraints 



We calculate a sample of mass spectra and check that this model provides a realistic mass 
spectrum consistent with current experimental bounds. 

Our numerical analysis is done in the following way. We fix the cutoff scale, M cut , which is 
of the same order as the KK scale, M 5 e~ kRn , from the relation (41) 

M - ~ * 3 x lol GeV = « * 1015 GeV ' 

We assume that contact term couplings between the MSSM fields and the SUSY breaking field 
in (42-50) are all 0(1) and adopt the following initial condition: 

M« /2 = 2 M x , (89) 
< u = m 2 Hd = M 2 X , (90) 

( m Q)ij = C Qij aiajM^ , 

(Q, a) -> (U,/3), (D, 7 ), (L,5), (E, e) , (91) 

A uij = -M x (Xn)ij + a uij PiajMx , (92) 

A dij = -M x {Y d ) ij + a dij jiajMx , (93) 

A eij = -Mx{Y e )ij + a eij e^jMx , (94) 

where M x was defined as Mx = | < F x > \/M 5 e~ kRlT and we set a natural range of the 
parameters as 0.1 < c^j,a^j < 1. The factor 2 in the right-hand side of (89) comes from the 
factor 4(g2) 2 in (42). Since M cut is around 10 15 GeV, SU(2) and SU(3) couplings of MSSM, 
9 , 4 =2 )fl , 4 =3 ) a ^ M cut take the value of 0.7. For simplicity, we fix the normalization of U(l) 
coupling at M cut as 0.7. Then we obtain the factor 2 in (89) from 

A{g2) 2 [fi r = M cut ] ~ 4-0.7 2 ~ 2 . 

Our aim is to prove that, in our 5D MSSM framework, there exists a mass spectrum that 
is consistent with the current experimental bounds. We arrange the parameters as 

1 for % J , 
c*ij, a*ij = 0.1 for i ^ j , 

to keep the flavor-violating terms as small as possible with a mild hierarchy among the model 
parameters. Now the model has three free parameters: 

Mx, 83 (~ #2 ~ 35i), tan/3 . 
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The KK scale is determined by £3 in (41). Since €3 is smaller than 0(1) from naturalness, (32) 
leads to the condition: 

1 £ <*3 £ m T /vcos(3 ~ 0.01/ cos (3. (95) 

Based on this setup, we calculate mass spectra for various values of (Mx, 83, tan/3) and 
check if they are consistent with the current experimental bounds, in particular, the lower 
bound of Higgs boson mass. With the flavor-violating soft terms predicted in our model, 
we then evaluate the rates of the lepton flavor violating processes, i.e. the branching ratios 
of /x — > e7 and r — > 7x7 decays, for each sparticle mass spectrum based on the technology 
developed in [20] and compared the results with current bounds. In our analysis of MSSM 
RG equations, we first ignore the off-diagonal terms in (91-94) from the initial condition and 
numerically solve the MSSM RG equations from M cut to low energies using Softsusy-3.1.4 [Zl] 
with Yukawa off-diagonal terms ignored. After this calculation, we add the off-diagonal terms 
to give the resultant spectrum. 

Below is the list of sample values of (Mx, 83, tan 0) that give realistic mass spectra consistent 
with the bounds on Higgs boson mass and \i — > ej and r — > /X7 branching ratios. We focus on 
the case with Mx < 600 GeV because light mass spectra are of more phenomeno logical interest. 
The lightest Higgs boson mass rrih, /x — > branching ratio (Br^), and r — > 7/7 branching ratio 
(Br T ) of each spectrum are also shown. The results shown here satisfy the current experimental 
bounds: m h > 114.4 GeV [22], Br^ < 1.2 x 10" 11 [23] and Br T < 4.5 x 10" 8 [21]. The entire 
mass spectra for three examples are shown in the Appendix. For larger values of £3 and/or 
tan/3, the spectrum violates the bound on /x — > branching ratio. For smaller Mx and/or 
tan/3, the lightest Higgs boson is too light. 
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To /x — > e7 process, loop diagrams containing the following terms contribute: 

< H° d > (A e ) 21 ~ e 2 5 1 v d M x ~ (5 1 /5 2 )m fl M x ~ ^m^M x , 
<H°> {A e ) 12 ~ e x 8 2 v d M x ~ (5 2 /5 l )m e M x ~ Sm e M x , 
K) 12 ~ 8 X 8 2 M 2 X ~ 1 -{53) 2 M X , 

(m|) 12 ~ e l e 2 M 2 x ~ • 
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The contributions from the terms (A e ) 21 , (A e ) 12 are almost independent of tan/3 and £3. For 
the mass spectra listed above, if (v4 e ) 21 , (A e ) 12 were the only source of lepton flavor violation, 
they would give Br^ ~ 7 - 8 x KT 12 for M x = 500 GeV and Br^ ~ 3 - 4 x 10~ 12 for M x = 600 
GeV. Hence we argue that, for the cases with small tan/3 and £3, the flavor-violating A-terms 
give dominant contributions. It is obvious that (A e ) 2 i contributes much more strongly than 
(A e )i2- On the other hand, the contribution from the term (m|)i 2 is sensitive to the values 
of tan/3 and £3, which is roughly proportional to (tan/3) 2 and (f^) 4 . (The net value of Br^ 
does not reflect this rule because of the interference between (A e ) 2 i contribution and {m 2 L )i2 
contribution.) We thus obtain the upper bounds on tan/3 and £3 when the contribution from 
the term (m|) 12 becomes dominant. The contribution from {m E )i2 is much suppressed by the 
tiny ratio m e /m T and has negligible impact on \x — > e'-f branching ratio. 
Here we summarize the features of the sample mass spectra listed above. 

(i) The typical SUSY breaking mass scale, Mx, can be as low as 500 GeV and the mass 
spectrum is within the reach of the LHC. 

(ii) The ratio £3/ tan/3 is around 0.01, which means that we need e 3 ~ 1 to have the tau 
Yukawa coupling. Therefore the 5D superfield of singlet tau is strongly localized towards the 
IR brane. 

(iii) \i — y branching ratio is always higher than O(10 -12 ) and the model can be tested by 
MEG experiment [25]. in the near future. 

(iv) t — > /i7 branching ratio is of the same order as [i — > cy branching ratio. 

The feature (ii) originates from the difference between the hierarchy of Si's and that of q's. 
The experimental bound on \i — > e'-f branching ratio constrains the terms (m|)i 2 and (m 2 E )±2, 
which are respectively proportional to <5i<5 2 and eie 2 , with the same extent. This gives a stronger 
limit on £3 than on €3 because <5j's have milder hierarchy. The orders of £3 and €3 are related 
through 

5 3 e 3 /tan/3 ~ m T /v ~ 0.01 , 

where tan (3 cannot be smaller than about 4 because otherwise the LEP II Higgs mass bound 
would not be satisfied. Therefore small £3 and large e% are favored in this model, which leads 
to the prediction that e% ~ 1. 

The feature (iii) results from the existence of the flavor- violating A-term, (A e )i2, and the 
fact that its contribution is independent of tan /3 and £3. The resultant branching ratio is within 
the future reach of MEG experiment. 

The feature (iv) is specific to this model because new physics models normally predict 
t — > fi^f branching ratio larger than [i — > e'-f branching ratio as new physics is more likely to 
affect 3rd generation particles than 1st and 2nd generations. This feature is a consequence 
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of the feature (ii). Since e 3 ~ 1, SU(2) singlet stau obtains large soft mass through gravity 
mediation on the IR brane and becomes a few times heavier than SU(2) singlet smuon and 
selectron if gravity mediation contributes positively as is usually assumed. For r — > fi'j process, 
the term (A e )s2 ~ e^5 2 Mx always contributes. One can compare its impact on Br T with that 
of (A e )2i on Br yi by comparing 

< H° d > {A e ) 32 I m T ~ M x vs. < H° d > (A e ) 21 / m„ ~ ^M x . 

The former is larger by the factor 3. However its effect is canceled by the larger mass of 
the singlet stau propagating in the loop diagram containing (As)32- The term (mf,^ also 
contributes when tan/3 is relatively large. (771^)32 is predicted to be about 3 times larger than 
(ml) i2- The contribution from (m 2 ) 32 mainly comes from two types of diagrams, SU(2) singlet 
smuon propagating in one diagram and singlet stau propagating in the other. However the 
latter is suppressed by the large stau mass, which partly cancels the effects of large (m|) 3 2. 
The term (m%) 2 3 contributes when tan /3 is relatively large. However, again, its contribution is 
suppressed by the large mass of the singlet stau propagating in the diagram. 

Finally, we discuss the prediction of our model on Am^ of K° — K° mixing and b — > 57 
branching ratio, based on the paper |26j . 

The following flavor- violating parameters predicted in our model are relevant to the K° — K° 
mixing: 
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Here we focus on the case with M x = 500 GeV. The average squark mass is around the same 
scale. Comparing the above predictions with Table 1 in [26], we see that, when tan /3 < 15, our 
predictions are below the limits that come from the experimental bound on Amx- 
For b — > 57 process, our model predicts the following flavor-violating parameters: 
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Again, we focus on the case with Mx = 500 GeV. Comparing the above predictions with Table 

6 in [26]. we find that, regardless of tan/3, our predictions are far below the limits that come 
from the experimental bound adopted by [26]. Even if we adopt the stronger bound in |27j . 
our predictions are still below the limits. 

7 Unusual NLSP and its flavor-violating decay 

In this section, we consider NLSP, which is mostly composed of SU(2) singlet charged sleptons. 
As we found in the previous section, our model favors e 3 ~ 1, i.e. the singlet stau superfield is 
localized towards the IR brane and receives large gravity mediation effects. This changes the 
flavor structure of charged slepton mass matrix. We here discuss the model's predictions on 
the flavor composition of NLSP. In the following, we fix 63 = 1. 
We first review the charged slepton mass matrix. Define 

A = A — ji tan (3 , 

where Ao indicates those parts of A-terms that are proportional to the corresponding Yukawa 
couplings. We denote RG contributions (gaugino mediation contributions) and D-term contri- 
butions to the soft masses of doublet selectron, smuon, stau and singlet selectron, smuon, stau 
by 

222222 

m Lli m L2i m L3i m Eli m E2i m E3 ■ 

Since we are focusing on NLSP, whose candidates are singlet sleptons, we neglect relatively 
small doublet slepton and doublet-singlet mixing terms. Further neglecting terms suppressed 
by m e /m r , we obtain the following approximate form of the charged slepton mass matrix: 
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(96) 

The factor c^ T denotes the coupling of the contact term among singlet smuon, singlet stau and 
the SUSY breaking sector, and is assumed to be 0(1). 

In the following, we consider two cases for completeness: 

(i) gravity mediation contributions to soft masses are positive, as is usually assumed. 
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(ii) they are negative. 

In case (i), SU(2) singlet stau becomes heavier than singlet smuon and selectron, so that 
the NLSP will be mainly composed of either singlet smuon or singlet selectron. In order to see 
which is the one, we consider the mixing mass between singlet smuon and singlet stau, namely 
(777,^)23. First note that the difference between m 2 El and m 2 E2 an d that between m L2 and m E2 , 
which arise from RG running, are given by 

16 , (M cut 



<i~ m E2 ~ i6^ ln h^J M ^ (^) tan 2 /3, (97) 
+ { --(cos 2 6 W - sin 2 6 w )Ml cos 2(3 + sin 2 6 W AfJ cos 2(3 } 



16tt 2 V. M w / 

where = m^/v and S* = rn 2 Hu — m? Hd + tr [wig — m^ — 2m 2 + m^ + m 2 ]. In deriving (98), we used 
the relation Mf, 2 ~ ^(g%) 2 M x at the scale M cnt shown in (42), and the fact that M 1 a /2 /(^) 2 is 
an RG invariant at the 1-loop level. 

Singlet smuon mixes with doublet smuon through the term Am^ and with singlet stau 
through the term c MT (m M /77i r )M|- . Solving these mixings, the mass eigenstate, which is still 
dominantly composed of singlet smuon, has mass estimated as 



m AM2». r 2 2(^m M ) 2 2(c /ir (m /i /m r )Mp 2 



cig — " b E2 T l ' 1V1 X 2 2 2 2 

m r m£ 2 - m E2 m E3 - m E2 



rn t \ n/r2 2 (/i tan/3 m M ) 2(c^ T {m^/m T )M 2 x ) 2 



E2 + 1 "~~ J -— — — ^-r-g^ ^ 



(24/16tt2) ln(M CTtt /M w ) ^M 2 - M 



A 



m| 2 + (1 - 2c 2 ) ( ^ ) M 2 X -mi tan 2 /3/^ . (99) 
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On the other hand, we can neglect the mixing terms between singlet selectron and other sleptons 
because they are suppressed by the tiny value of m e /m T . 

We now compare the mass of singlet selectron m 2 El with the mass of the eigenstate, m 2 ig . 
Using (97) and (99), we have 

m% -m 2 E1 * In M 2 X [y,f tan 2 /3 + (1 - 2c 2 r ) (^j M 2 X - m 2 tan 2 fi/fi 



Til 



m, 



2 



l ~Kr) [-^ M x . (100) 



In summary, in the determination of the NLSP mass, the effect of the mixing term between the 
singlet smuon and singlet stau dominates over the Yukawa RG contributions and the effects of 



27 



singlet-doublet mixing terms. NLSP is singlet-smuon-like for c^ T > 1/ y/2, whereas it is singlet- 
selectron-like for a relatively small coupling c Mr < l/y/2. Since the factor c^ T affects r — > fi'j 
branching ratio through the term (m E )2z = c tlT {m il /m T )M'j c , we have a connection between the 
value of Br T and the flavor of NLSP. 

If NLSP is smuon-like, it decays mainly into \x and gravitino. However our model predicts 
that NLSP can contain a considerable amount of stau component due to the mixing term 
(m 2 E ) 2 3, which can be as large as [m il jm 1 )M\. Therefore the branching ratio of a flavor- 
violating NLSP decay into r and gravitino can be as large as 

Br(NLSP r ^3/2) ~ K/m r ) 2 ~ ^ . (101) 

On the other hand, in the context of MFV, the mixing terms of singlet smuon and other sleptons 
are much smaller as we evaluated in section 5, and such flavor- violating NLSP decays are much 
suppressed. Thus the flavor- violating decay of NLSP provides a distinct signature of our model. 
If NLSP is selectron-like, it decays mostly into electron and gravitino. The model predicts that 
the mixing terms of singlet selectron and other sleptons are suppressed by m e , although they 
are much larger than in other models with MFV. Still, as a distinct signature of our model, we 
expect to observe a rare NLSP decay into r and gravitino with the branching ratio as large as 

Br(NLSP ->■ r tp 3/2 ) ~ (m e /m r ) 2 ~ 1(T 7 . (102) 

Note that the lifetime of NLSP is estimated as 

I < Fy > I 2 MlM 2 
t NLSP ~ 48vrJ ^ ~ 48tt . x c % , 

m NLSP {m N LSPr 

which is ~ 10- 3 sec for m NLS p = 300 GeV, M x = 500 GeV and M cut = 10 13 GeV (5 3 ~ 0.1). 
Such a long-lived NLSP, once produced at the LHC, decays outside the detector. There have 
been interesting proposals [28] for the way to trap charged NLSPs outside the detector. Detailed 
studies of the NLSP decay can allow us not only to measure the gravitino mass and the four- 
dimensional Planck mass but also to test the flavor-structure of NLSP predicted in our model. 

In case (ii), the gravity mediated contributions are negative, so that SU(2) singlet stau is 
the dominant NLSP component. Although this singlet-stau-like NLSP is as usual in SUSY 
models with gravitino LSP, we again expect to observe a rare decay of NLSP as a signature of 
the model. The branching ratio of NLSP decay into \x and gravitino is predicted to be as large 
as 

Br(NLSP -> fi ^3/2) ~ (mJm T f ~ -L . (103) 
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8 Conclusion 



We have investigated a simple 5D extension of MSSM in RS spacetime, where 5D geometry 
controls both the SUSY breaking mediation mechanism and the Yukawa coupling hierarchy. 
The Yukawa coupling hierarchy is naturally explained by the localization of matter superfields 
in the 5D bulk. SUSY breaking effects arise from two sources: contact terms between the SUSY 
breaking sector and the MSSM fields (gravity mediation), and RG effects (gaugino mediation). 
The former are flavor-violating and the latter flavor-conserving. Using the experimental data 
on fermion masses and mixings, we have determined the 5D disposition of matter superfields 
and calculated SUSY breaking mass spectra including flavor-violating terms. We have numer- 
ically checked that our framework can give a realistic mass spectrum consistent with all the 
experimental constraints. 

We have estimated the flavor- violating effects induced by RG running in the context of MFV 
and compared them with those from gravity mediation predicted by the model. We have proved 
that our model provides a different pattern of flavor-violating terms, namely, flavor violation 
of A-terms and SU(2) singlet soft masses can be much larger than the MFV case. 

Our model has several distinct predictions. First, the \i — > branching ratio is larger 
than O(10 -12 ), regardless of tan/3 and the seesaw scale, for sparticle masses < 2 TeV. This 
originates from the basic structure of our model, namely, the hierarchy of Yukawa couplings 
and the gravity mediation contributions to A-terms are rooted on the same 5D disposition of 
the matter superfields. Hence, using the experimental data on the charged SM fermion masses, 
CKM matrix and the neutrino oscillation parameters, we can fix the orders of flavor-violating 
A-terms. Second, r — > /X7 branching ratio may not be larger than /x — > c-f branching ratio. This 
is because our model predicts IR-localized singlet stau superfield and it gains an additional soft 
mass through gravity mediation on the IR brane. Third, since RS geometry warps down the 
effective cutoff scale, which can be characterized by the seesaw scale, gravitino is LSP and the 
dark matter candidate. Forth, NLSP is likely to be either smuon-like or selectron-like because 
the gravity mediation contribution pushes up the singlet stau mass. Furthermore our model 
predicts flavor- violating NLSP decays with the rates much higher than those expected in the 
MFV case. 

The Yukawa coupling hierarchy is one of the long-standing problems in the Standard Model. 
The 5D MSSM on Randall- Sundrum background offers a solution to this problem from the 
geometrical point of view. In this model, flavor-violating soft SUSY breaking terms have the 
same geometrical origin as the Yukawa coupling hierarchy. Therefore, the flavor structure 
among sparticles can be a clue to understand the origin of flavors among the SM fermions, even 
if the origin lies at an energy scale far above the electroweak scale. Gravitino is always LSP due 
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to the warped geometry, and NLSP is predicted to be dominantly composed of SU(2) singlet 
sleptons and long-lived. At collider experiments, in this case, supersymmetric events can be 
fully reconstructed without missing energy, which allows us not only to identify the dominant 
flavor content of NLSP, but also to measure the rates of flavor-violating decays of sparticles if 
it is sizable. Thus, our framework can be tested in the future. 
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Appendix: Sample Mass Spectra 
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Table 1: Particle mass spectra for three samples for different values of (Mx, tan/3, 5). 
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